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Macroscopic Observability of Spinorial Sign Changes: A Reply to Gill 

Joy ChristiarQ 

Einstein Centre for Local-Realistic Physics, 15 Thackley End, Oxford 0X2 6LB, United Kingdom 

In a recent paper Richard Gill has criticized an experimental proposal published in a journal of 
physics which describes how to detect a macroscopic signature of spinorial sign changes. Here I point 
out that Gill’s worries stem from his own elementary algebraic and conceptual mistakes, and present 
several event-by-event numerical simulations which bring out his mistakes by explicit computations. 


In a recent paper a mechanical experiment has been proposed to test possible macroscopic observability of spinorial 
sign changes under 2n rotations (TJ. The proposed experiment is a variant of the local model for the spin-1/2 particles 
considered by Bell 2], which was later further developed by Peres providing pedagogical details |0. Our experiment 
differs, however, from the one considered by Bell and Peres in one important respect. It involves measurements of the 
actual spin angular momenta of two fragments of an exploding bomb rather than their normalized spin values, ±1. 

It is well known that angular momenta are best described, not by ordinary polar vectors, but by pseudo-vectors, or 
bivectors, that change sign upon reflection 01- One only has to compare a spinning object, like a barber’s pole, with 
its image in a mirror to appreciate this elementary fact. The mirror image of a polar vector representing the spinning 
object is not the polar vector that represents the mirror image of the spinning object. In fact it is the negative of 
the polar vector that does the job. Therefore the spin angular momenta L(a, A) in the theoretical analysis of the 
experiment proposed in the paper jl] have been represented by bivectors using the powerful language of geometric 
algebra @. They can be expressed in terms of the bivector basis (which are graded basis) satisfying the sub-algebra 

L P (X) L U (X) = — g P i/ — 'y ' e pvp L P {X), (1) 

p 

where A = ± 1 represents a choice of orientation of a unit 3-sphere T|. This brings us to the first of several elementary 
mistakes made by Gill. In his paper [3] he claims that, from the above equation, using L p {X = —1) = — L p {X = +1), 

Lp{X = +1) L u {X = +1) — Lp{X = —1) T„(A = —1) = — Cp Up L p {X = +1) + e pvp L p (X = — 1) (2) 

p p 

implies 0 = -2^ e pvp L p {X = +1) = +2^ e pup L p (X = -1), (3) 

p p 

which in turn implies L p {X = +1) = L p {X = —1) = 0. (4) 

It is not difficult to see, however, that this claim is manifestly false. Gill’s mistake here is to miss the summation over 
the index p. And since the basis elements L p {A = +1) and L p (X = —1) do not vanish in general, contrary to his claim 
neither the spin angular momentum L(a, A) about the direction a nor the spin angular momentum L(b, A) about the 
direction b vanish in general. Only the spin angular momentum L(a x b, A) about the mutually orthogonal direction 
a x b vanish, as in Eq. (110) of the paper pj. This is of course consistent with the physical fact that there is no third 
fragment of the bomb spinning about the direction a x b exclusive (as well as orthogonal) to both directions a and 
b. To see that only the spin angular momentum L(a x b, A) about the direction a x b vanishes, all one has to do is 

to contract Eq. m above with the vector components a p and b v , on both sides, and then follow through the steps in 

Eqs. ([2]) and ([3]). It is also crucial to appreciate that the spin angular momenta L(s, A) (ie., the bivectors) trace out 
an su(2) 2-sphere within the group manifold SU(2) ~ S 3 , not a round S 2 within IR 3 as Gill has incorrectly assumed. 

Unfortunately the conceptual mistakes made by Gill are even more serious than the above elementary mathematical 
mistakes [7] @. As already noted, what are supposed to be measured in the proposed experiment are the spin angular 
momenta L(a, A) and L(b, A) themselves, not their normalized spin values ±1 about some directions a and b, where 

L(a, A) s A(a, A) = A(e x A e y A e z ) • a = ±1 spin about a (5) 

and L(b, A) s B(b, A) = A(e x A e y A e z ) • b = ±1 spin about b. (6) 

This is in sharp contrast to what are measured as dynamical variables in the model considered by Bell and Peres [i[. 


Electronic address: jjc@alum.bu.edu 



2 


In practice the above dynamical variables are supposed to be measured in the proposed experiment by directly 
observing the polar vectors s fc dual to the random bivectors L(s, A fc ). It is only after all the runs of the experiment are 
completed and the vectors s fc are fully recorded, the traditional dynamical variables sign ( + s fc • a) and sign (— s fc ■ b) 
supposed by Bell are to be calculated, by an algorithm extraneous to the actual experiment. These calculations may 
be done, for example, years after the experiment has been completed. Thus they are not supposed to be an integral 
part of the physical experiment itself, and this fact must be accounted for in the ensuing statistical analysis of data. 

Since Gill has failed to appreciate the difference between the actual dynamical variables in the proposed experiment, 
namely the bivectors L(a, A fc ) and L(b, X k ), and the subsequently calculated traditional variables sign (+ s k • a) and 
sign (— s k ■ b), he has ended up misapplying his calculation, which leads him to a result not relevant to the problem at 
hand. In order to correctly calculate the correlation between the variables sign ( + s fc • a) and sign ( — s k ■ b) one must 
first specify their relation to the actual dynamical variables L(a, A fc ) and L(b, A fc ) to be observed in the experiment: 

S' 3 9 ±1 = sign( + s k • a) = gf(a, X k ) = lim { - D(a) L(s, A fc )} = J + ] ^ \ + * (7) 

s —a I — 1 II A = — 1 

and 

S 3 9Tl = sign(-s k -b) = 38(b, X k ) = lim {+D(b)L(s, A fe )} = * u =+ ! (8) 

s-tb 1+1 II A = — 1 , 

where the orientation A of S 3 is assumed to be a random variable with 50/50 chance of being +1 or — 1 at the moment 
of the bomb-explosion, making the spinning bivector L(a, A) a random variable relative to the detector bivector D(a): 

L(a, A) = { 0#1 L M (A) } = A{ a v D v } ee A D(a). (9) 


From the above discussion it should be clear that the raw scores sign (+ s k ■ a) and sign (— s k ■ b) would be 
generated in the experiment with different bivectorial scales of dispersion, or different standard deviations, as explained 
between Eqs. (102) and (109) in Ref. QJ. Therefore the calculation of the correlation between these raw scores, as well 
as the derivation of the Tsirel’son’s bound on the strength of possible correlations, must be carried out with some 
care. In fact there are at least four physical considerations one must be mindful of before proceeding further: (1) the 
spin angular momenta are represented, not by polar vectors, but by bivectors that change sign upon reflection; (2) 
scalars and bivectors, despite being elements of different grades, are treated on equal footing in geometric algebra; (3) 
what would be actually observed in the proposed experiment are not the raw scores ^/(a, X k ) = sign ( + s k ■ a) and 
38 (h, X k ) = sign ( — s k • b) but the standard scores A{ a, A fc ) = L(a, X k ) and B{ b, A fc ) = L(b, A fc ); and (4) the correct 
association between the raw scores (a, X k ) and 38{h, X k ) can be inferred only by calculating the covariance between 
the corresponding standard scores A(a, X k ) and B( b, X k ). With these in mind, let us consider four reference vectors a, 
a', b, and b'. Then the bound on the corresponding CHSH string of expectation values [E}, namely, on the coefficient 

f (a, b) + £(a, b') + S(a!, b) - £{a!, b'), (10) 

can be derived using the four joint expectation values of the raw scores — such as srf (a, A) and 38 { b, A) — defined as 


£(a, b) 


lim 

n»l 


1 

- V s/(sl, X k ) 381 b, A fe ) 
n 

fe = l 


(^a(A)^b(A)). 


( 11 ) 


This allows us to express the above CHSH string of expectation values simply as a string of four averages as follows: 

(*4(A)# b (A)) + (^ a (A)^b'(A)) + (^a'(A)^b(A)} - (^(A)^b'(A)}- (12) 

It is in the next step that Gill makes his gravest mistake. He surreptitiously replaces the above string of four separate 
averages of numbers that are generated with different bivectorial scales of dispersion with the following single average: 

(*4(A)# b (A) + <(A) J? b '(A) + *4'(A)# b (A) - *4/(A)# b '(A)). (13) 

As innocuous as this step may seem, it is in fact an illegitimate mathematical step within the context of the experiment 
proposed in Ref. [ij. But this illegitimate maneuver does allow Gill to reduce the above average at once to the average 

^ a (X){38 b (X) + 38 w (X)} + tfk'(A){#b(A)-#b'(A)}). 


(14) 
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And since ^b(A) = ±1, if |^b(A) + ^b'(A)| = 2, then |^b(A) — &b'(X)j = 0, and vice versa. Consequently, using 
£/ a (X) = ±1, it is easy to conclude that the absolute value of the above average cannot exceed 2, as Gill has concluded. 

As compelling as this conclusion by Gill may seem at first sight, it is entirely false. It is based on his illegitimate and 
careless maneuver of replacing the string of four separate averages of random variables m with a single average 0. 
Such a move is justified on mathematical grounds only if all four random variables g/ a (X), & b (X), ^'(A), and <^b'(A) 
are on equal statistical and geometrical footings. But as we noted earlier, each one of these variables is generated 
with a different bivectorial scale of dispersion, or a different standard deviation, as explained between equations (102) 
and (109) in Ref. [1]. Therefore the theoretical prediction of the correlation m between these raw scores, as well as 
the theoretical derivation of the Tsirel’son’s bound on the strength of possible correlations, must proceed as follows. 

Since the correct association between the raw scores stf (a, A) and &(b, A) can be inferred only by calculating the 
covariance between the corresponding standardized variables A( a, A) and B( b, A) as calculated in Ref. Q, namely 

A a (A) = A(a, A) = L(a, A) (15) 

and B h (X) = B( b, A) = L(b, A), (16) 

the correlation between the raw scores srf (a, A) and A) must be obtained by evaluating their product moment 


£(a, b) 


lim 

1 


1 n 

- J2 AL(a, X k ) B( b, A fc ) 


(17) 


The numerical value of this coefficient is then necessarily equal to the value of the correlation calculated by Eq. GD- 

Il sing the above expression for £ (a, b) the string of expectation values (flOl) can now be rewritten as a single average 
in terms of the standard scores A a (A) and l?b(A), because now they are on equal statistical and geometrical footings: 


lim 

n »1 


1 

- V {AffX k ) B b (X k ) + A a (X k )B h ,(X k ) + A a /(A fe ) B h (X k ) 


A a ,(X k )B w (X k )} 


(18) 


But since A a (A) = L(a, A) and B b (X) = L(b, A) are two independent equatorial points of S 3 , we can take them to 
belong to two disconnected “sections” of S 3 [ie., two disconnected su(2) 2-spheres within S 3 ~ SU(2)], satisfying 

[A n (A), B n ' (A)] = 0 V n and n' £ IR 3 , (19) 

which is equivalent to anticipating a null outcome along the direction n x n' exclusive to both n and n'. If we now 
square the integrand of equation (11811 . use the above commutation relations, and use the fact that all unit bivectors 
square to —1, then the absolute value of the Bell-CHSH string (fllil) leads to the following variance inequality [j|: 


|f (a, b) + £(a, b') + £(a!, b) 


£(a', b')| 




lim 

n » 1 


1 

- V {4 + 4^ aa /(A fe ) ^b' b (A fc ) } 

[ n fc = i 


( 20 ) 


where the classical commutators 


f^aa'(A) — [A a (A), A a /(A)] — — A axa '(A) (21) 

and 

^b'b(A) := i [Bb-(A), B b (A)] = — -Bb'xb(A) (22) 

are the geometric measures of the torsion within S 3 [EJ. Thus, as discussed in the paper, it is the non-vanishing torsion 
X7 within the 3-sphere—the parallelizing torsion which makes its Riemann curvature vanish—that is responsible for 
the stronger-than-linear correlation. We can see this from Eq. (EU1) by setting £7 = 0, and in more detail as follows. 

Using definitions m and (fl6l) for A a (A) and B b (X) and making a repeated use of the well known bivector identity 

L(a, A) L(a', A) = - a-a' - 


L(a x a', A) 


(23) 
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the above inequality can be further simplified to 


|£(a, b) + £(a, b') + £(a', b) - £(a', b')| < 



4 — 4 (a x a') • (b' x b) — 4 lim 

n »i 

^ 2 \/l-(axa')-(b'xb) — 0 , 


D(z) 


k = 1 


(24) 


where z = (a x a') x (b' x b). The last two steps follow from the relation ([9|) between L(z, A) and D(z) and the fact 
that the orientation A of S 3 is evenly distributed between +1 and —1. Finally, by noticing that trigonometry dictates 


— 1 < (a x a') • (b' x b) < +1. 


(25) 


the above inequality can be reduced to the form 

| £(a, b) + £(a, b') + £(a', b) - £(a', b') | ^ 2y/2, (26) 

exhibiting the correct upper bound on the strength of possible correlations. Thus the stronger bounds of —2 and +2 
calculated naively by Gill in his paper are simply incorrect. More importantly, as correctly done in the paper [l|, the 
correlation function for the bomb fragments respecting the SU(2) Lie algebra su(2) can indeed be calculated to yield 


£(a, b) = lim 

n » 1 


1 " 

— ^ {sign (+s fc • a)} {sign (—s fc • b)} 
. fc = l 


— a • b, 


(27) 


where n is the total number of trials performed. Therefore the remaining comments by Gill are anything but justified. 

In fact there is considerable confusion in Gill’s attempted misrepresentation of the proposed experiment [7j. Rather 
surprisingly, his comments fail to distinguish between the theoretical derivation of the correlation presented in Eq. (110) 
of Ref. [l] and the practical calculation of the correlation expressed in the above equation in terms of the raw scores. For 
example, he asserts that “correlations should only be computed the usual way using actual experimental outcomes...” 
This reveals that Gill has either not bothered to read the experimental proposal discussed in paper [I), or has failed 
to understand the crucial difference between the traditional Bell-type experiments and the one proposed in the paper. 

The important question here is: What would be the “actual experimental outcomes” in the proposed experiment? 
From the above discussion and the discussion in the Section 4 of Ref. flj it is quite clear that what would be actually 
observed in the experiment are the bivectors L(a, A) and L(b, A) about the directions a and b, respectively. The 
raw scores sign (+ s fc ■ a) and sign ( — s k ■ b) may then be calculated, but by an algorithm extraneous to the actual 
experiment, and long after the experiment has actually been completed. Thus they would not be an integral dynamical 
part of the physical experiment itself. On the other hand, thanks to the precise geometrical relation between the raw 
scores sign ( + s k • a) and sign ( — • b) and the standard scores L(a, A) and L(b, A) given by the equations © and 

©, the statistical association between the raw scores can nevertheless be inferred by calculating the covariance of the 
corresponding standardized variables A a (A) and i?b(A) as demonstrated above, thus predicting the correlation (1271) . 

Finally, Gill falsely and disingenuously claims that Section 5 of Ref. [lj reproduces a sign error which is also present 
in my earlier work Q. There is in fact no such error, as explained already in Ref. Q and in several chapters of Ref. [jf. 

In conclusion, the criticism of the proposed experiment by Gill is a travesty. No physicist should be deceived by it. 


Appendix A: Refutations of Gill’s Mistaken Claims by Explicit Numerical Simulations 


The above refutations of Gill’s fallacious claims have been independently verified by Albert Jan Wonnink [9j in an 
explicit numerical simulation of Eq. m, by means of a specialized program for geometric algebra based computations. 
In other words, Gill’s mistaken claims have been refuted by Wonnink by numerically computing the expectation value 


£(a, b) 


lim 

n » 1 


n 


E 


L(a, A fe ) L(b, A fc ) 


n 


= —a • b. 


(Al) 
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To understand this computation, recall that L(a, X k = +1) = +1 ■ a and L(b, X k = —1) = — I ■ b represent the two 
spins of the bomb fragments, where I := e x A e y A e z with +1 representing the right-handed orientation of S 3 and —I 
representing the left-handed orientation of S 3 . Consequently we may consider the following two geometric products, 

( + I-a)( + /-b) = -a-b - ( + /)• (ax b) (A2) 

and 

(— J• a)( — J• b) = - a-b - (-/)• (ax b), (A3) 

as discussed in my earlier reply to Gill (cf. Eqs. (15) and (16) of Ref.@). The two possible orientations of S 3 may 
then be thought of as the random hidden variables A = ± 1 (or the initial states A = ± 1) of the two bomb fragments. 

Now in traditional geometric algebra as well as in the GAViewer program 0 employed by Albert Jan Wonnink the 
volume form of the physical space is fixed a priori to be +1, by convention. This convention is inconsistent with the 
physical process of spin detections delineated in the Eqs. 0 and 0 above. Therefore a translation of the geometric 
product L(a, X k = —1) L(b, X k = —1) = (— I ■ a)( — I ■ b) for the GAViewer built on the right-handed form +1 is 
necessary, which can be inferred from Eqs. (IA2I) and (IA3I) by recalling that bxa = — axb and rewriting Eq. (IA2I) as 

(+ / • b)( + / • a) = — a ■ b + ( + /) • (a x b) = —a-b — ( — J)-(axb). (A4) 

Comparing the right-hand sides of Eqs. (1A3I) and (1X41) it is now quite easy to recognize that the desired translation is 

( — / • a)( — J • b) —> (+ / • b)( + / • a). (A5) 

In terms of this translation the geometric products appearing in the expectation function m can be expressed as 

L(a, X k = +1) L(b, X k = +1) = ( + / • a)( +1 ■ b) (A6) 

and L(a, A fc =-1) L(b, A fc =-1) = (+ / ■ b)( + / • a). (A7) 

In other words, when X k happens to be equal to +1, L(a, A fc ) L(b, A fc ) = (+1 ■ a)( +1 ■ b), and when X k happens to 
be equal to —1, L(a, X k ) L(b, X k ) = ( +1 ■ b)( +1 ■ a). Consequently, the expectation value (IA1I) reduces at once to 

£(a, b) = ^ (+1 ■ a) (+1 ■ b) + ^ (+1 ■ b) (+1 ■ a) = — a • b + 0, (A8) 

because the orientation A of S 3 is necessarily a fair coin. Here the last equality follows from the Eqs. (IA2D and (1A4D . 
Given the translation (1A5I) , it is now easy to understand the simulation code of Ref. ;9j , with its essential lines being 

if A = +1, then add ( +1 ■ a)( +1 ■ b), (A9) 

but if A = —1, then add ( + / ■ b)( +1 ■ a). (A10) 

Not surprisingly, the expectation value computed in this event-by-event simulation prints out to be £(a, b) = — a • b. 

In complement to the above unambiguous demonstration, it is also possible to verify the correlation (IA8I) numerically 
in an event-by-event simulation within a rum-Clifford algebraic representation of S 3 , as done, for example, in Ref. 11]. 

Thus, once again, the criticism of the proposed experiment by Gill is exposed to be vacuous by explicit computations. 


Appendix B: Correlations Among the Raw Scores is Equal to the Covariance Among the Standard Scores 


For the completeness of the arguments presented above, in this appendix let us prove the following equality explicitly: 


£(a, b) = lim 

n » 1 


It 

- i si 9 n (+s fc • a )} {sign- (-s fc • b )} 


fc = l 


= lim 

i 


J2 L(a, A fe ) L(b, A fc ) 


k = 1 


= —a • b. 


(Bl) 


We begin with the central equalities in the prescription of the remote measurement events given in Eqs. 0 and 

sign( + s k ■ a) =&/(&, X k ) = lim { — D(a) L(s, X k ) } (B2) 

and sign( — s k -b) = J?(b, X k ) = lim { + D(b) L(s, X k ) } , (B3) 

s-tb 1 ’ 
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where the same scalar number si (a, X k ) = ±1 is expressed in two different ways — as a grade-0 number sign ( + s k • a) 
on the left hand, and as a product — D(a) L(a, A fc ) of two grade-2 numbers, — D(a) and L(a, A fc ), on the right hand. 

Now the correlation between si (a, A fc ) and &(b, \ k ) can be quantified by the product-moment correlation coefficient 


£(a, b) 


lim 

n> 1 


1 

- E {*>> A ") - xk ) - 

. fc = 1 _. 

a (si) a(38) 


(B4) 


where si( a, A) and 38(b, A) are the average values of si and 38 and (7 (si) and a (38) are their standard deviations: 


a (si) 


1 " 1 

ly 

Tl ^' 1 


2 

1 ” i 


si (a, X k ) — si( a, A) 

and a (38) = 

- E 

n ^ J 1 

38( b, X k ) - 38( b, A) 

fc=l 


\ 

k = 1 



(B5) 


Accordingly, let us first consider si (a, X k ) — sign ( + s k - a) and 38( b, A fc ) = sign (— s k ■ b) from the left sides of the 
equalities (1B2I1 and (11 hit) . Written in this form, these variables cannot be factorized into products of a random number 
and a non-random number. We are therefore forced to treat them as irreducible random variables. Moreover, since A 
is a fair coin, s has equal chance of being parallel and anti-parallel to a. Consequently, it is easy to work out from the 
above formulae that in the present case si( a, A) = 38( b, A) = 0 and a (si) = a(38) = 1, which immediately gives us 


£(a, b) = lim 

n» 1 


1 

- E { si 9 n (+ sfe ' a )} i si 9 n (-s fe ’ b )} 


(B6) 


Thus, we recognize that the above expression of the observed correlations — which is traditionally employed by the 
experimentalists — is simply a special case of the Pearson’s product-moment correlation coefficient defined in Eq. (IB4I) . 

Next, consider si (a, X k ) = — D(a) L(a, X k ) and 38( b, A fc ) = + D(b) L(b, X k ) from the right sides of the equalities 
(IB2I) and (IB3I) . Written in this form, these very same variables are now factorized into geometric products of a random 
number and a non-random number. Evidently, the variable si( a, X k ) is a product of a random spin bivector L(a, X k ) 
(which is a function of the random variable A) and a non-random detector bivector — D(a) (which is independent of the 
random variable A). In other words, the randomness within si( a, X k ) originates entirely from the randomness of the 
spin bivector L(a, X k ), with the detector bivector — D(a) merely specifying the scale of dispersion within L(a, X k ). 
Consequently, as random variables, s/( a, X k ) and 38( b, A fc ) are generated with different standard deviations, or 
different sizes of the typical error. s/( a, X k ) is generated with a typical error — D(a), whereas 38( b, X k ) is generated 
with a typical error +D(b). And since errors in linear relations propagate linearly, the standard deviation a(si) of 
si (a, X k ) is equal to — D(a) times the standard deviation of L(a, A fe ), and the standard deviation a(38 ) of 38 ( b, X k ) 
is equal to + D(b) times the standard deviation of L(b, A fc ), which can be worked out using formulae similar to (IB5I) : 


& (si) = — D(a) er{L(a, A fc )} 
and a(38) = +D(b) cr{L(b, A fc )} 


-D(a) 

+ D(b) 


\ n 

\ k = l 


E I L(a, X k ) - L(a, A) 


\ n ' 
\ k = l 


E I L(b, X k ) - L(b, A) 


— D(a) 

+ D(b). 


(B7) 

(B8) 


Here I have used 11 L(a, X k ) 11 = 11 L(b, X k ) 11 = 1 because all L(s, A fc ) are unit bivectors, and L(a, A) = L(b, A) = 0 
on the account of A being a fair coin. These standard deviations are derived much more rigorously in Refs, [l] and Q. 

If we now substitute si (a, X k ) = — D(a) L(a, X k ) and 38 ( b, A fc ) = + D(b) L(b, X k ) into the same formula (IB4I) used 
previously for the product-moment correlation coefficient Q, together with si (a, X) = 38 ( b, A) = 0, a(si) = — D(a), 
a(38) = +D(b), {si( a, X k ) / a(si)} = L(a, A fc ), and {38 ( b, A k )/a(38)} = L(b, A fc ), then it immediately reduces to 


£(a, b) = lim 

1 


E L(a, A fc ) L(b, X k 


= 1 


(B9) 


Consequently, putting the results 


and 


together, we finally arrive at the equality we set out to prove: 


£(a, b) = lim 

n> 1 


1 

- E { si 9 n (+s fc ■ a )} {sign (-s fc • b)} 

_k = 1 


= lim 

n» 1 


1 " 

- E L( a > Afe ) b(b, X k ) 

n z ' 


—a ■ b. 


(BIO) 
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Appendix C: A Simplified Local-Realistic Derivation of the EPR-Bohm Correlation 12 


As in Eqs. (JTJ) and let the spin bivectors q= L(s, A fc ) be detected by the detector bivectors D(a) and D(b), giving 
S 3 9 trf (a, X k ) := hm (-D(a)L(s, A fe )} = j + J tf A* = - 1 } ^ ^ ) = ° (C1) 


and 


- 1 if X k = + 1 


S 3 9 #(b, A fc ) := hm {+L(s, A fe )D(b)} = < + l i{ X k = _ 1 \ with \^( b > A ") / = 0, (C2) 


where the orientation A of S 3 is assumed to be a random variable with 50/50 chance of being +1 or — 1 at the moment 
of the pair-creation, making the spinning bivector L(n, X k ) a random variable relative to the detector bivector D(n): 

L(n, A fc ) = A fc D(n) «<=► D(n) = A fe L(n, A fc ). (C3) 

The expectation value of simultaneous outcomes ,c/(a, X k ) = ±1 and X k ) = ±1 in S 3 then works out as follows: 


£(a, b) = lim 


1 " 

- V s/Isl, X k ) &(b, X k ) 

n ' 


= lim 

n —y oo 


= lim 


= lim 

n —> oo 


= lim 

n —> oo 


= lim 


k = 1 


E 

*:=i 


within S 3 := the set of all unit (left-handed) quaternions [l([5j (C4) 


lim { -D(a)L(s, A fe )} 


lim { + L(s, A fc ) D(b)} 


s —y b 


- E Jim {-D(a)L( S , X k ) L(s, A fc )D(b) = q(a, b; s, A fc )} 


k'= 1 s-tb 


(conserving total spin = 0) (C5) 
(by a property of limits) (C6) 


s lim { - A fc L(a, A fc ) L(s, A fc ) L(s, X k ) A fc L(b, X k )} 


fc = 1 s-^b 


(all bivectors in the spin basis) (C7) 


s lim {-L(a, X k ) L(s, X k ) L(s, X k ) L(b, A fe )} 


k — 1 s — y b 


(scalars X k commute with bivectors) (C8) 


1 " 

-E L ( a ’ A ") L ( b ’ A ") 

71 ' 


k = 1 


(follows from the conservation of zero spin angular momentum) (C9) 


= a • b hm 


= a ■ b — hm 

n —y oo 


E L ( axb > Afc ) 


k = 1 


1 n 

~y x k 

k = 1 


(NB: there is no “third” spin about the direction a x b) (CIO) 


D (a x b) (summing over counterfactual detections of “third” spins) (Cll) 


= — a • b + 0 (because the scalar coefficient of the bivector D(a x b) vanishes in the n —> oo limit) (C12) 

Here the integrand of (C6) is necessarily a unit quaternion q(a, b; s, A fe ) since S 3 remains closed under multiplication; 
(C7) follows from using (C3); (C8) follows from using A 2 = +1; (C9) follows from the fact that all unit bivectors such 
as L(s, X k ) square to —1; (CIO) follows from the geometric product K3l) : (Cll) follows from using (C3); (C12) follows 
from the fact that initial orientation A of S 3 is a fair coin; and (C6) follows from (C5) as a special case of the identity 


lim { — D(a) L(s, A fc )| lim { + L(s, A fe ) D(b)| 

— >a' K J s^h' K } 


s —> a 
s —y 


= lim - D(a)L(s, A fc )L(s, A fe )D(b) 


(C13) 


which can be easily verified either by immediate inspection or by recalling the elementary properties of limits. Note that 
apart from the assumption (C3) of initial state A the only other assumption needed in this derivation is the conservation 
of zero spin angular momentum. These two assumptions are necessary and sufficient to dictate the singlet correlation: 


£(a, b) = hm 


1 


E ^( a > A ") #( b > Xk ) 


k = 1 


= — a ■ b . 


(C14) 


This demonstrates that EPR-Bohm correlations are correlations among the scalar points of a quaternionic 3-sphere. 




























Appendix D: On the Fatal Mistake Made by John S. Bell in the Proof of His Famous Theorem 


Consider - in slightly modern terms - the standard EPR-Bohm type experiment envisaged by John S. Bell to prove 


his famous theorem . Alice is free to choose a detector direction a or a' and Bob is free to choose a detector direction 
b or b' to detect spins of the fermions they receive from a common source, at a space-like distance from each other. 
The objects of interest then are the bounds on the sum of possible averages put together in the manner of CHSH [l3 |. 


£(a, b) + £(a, b') + £{ a', b) - £(a', b'). 


with each average defined as 


£(a, b) = lim 

n> 1 


- £ ^(a, A fc ) mh, X k ) 

n ^ 


n , 

L k= 1 


= (<*4(a)^fc(b) 


(Dl) 


(D2) 


where ,c/(a, X k ) = ^4(a) = ±1 and J?(b, X k ) = & k (b) = ±1 are the respective measurement results of Alice and Bob. 
Now, since &/ k (a) = ±1 and B§ k {b) = ±1, the average of their product is —1 < ^ ^4(a) 3§ k { b) ^ ^ +1. As a result, 
we can immediately read off the upper and lower bounds on the string of the four averages considered above in CEO): 


- 4 < 


■s4(a)^ fc (b) ^ ^ ^4(a) 3§k(W) ^ + (&/ k (a!) @ k (b)^ - ^ M:(a') ^fe(b') ^ < +4. (D3) 


This should have been the final conclusion by Bell. But he continued. And in doing so he made one of his gravest 
mistakes. He replaced the above string of four separate averages of binary numbers with the following single average: 

£(a, b) + £(a, b') + £(a', b) - £(a', b') —■» ^ ,*4 (a) ^ (b) + <a4(a)^ fc (b') + &Z k (a.') ^ fc (b) - ^4(a') ^fc(b') ^ . 

(D4) 

As innocuous as this step may seem, it is in fact an illegitimate step, for what is being averaged on the RHS are 
un-observable, and hence un-physical quantities. But this illegitimate step allows one to reduce the above average to 

(^fc(a) { # fc (b) + # fc (b') } + *St(a') {@ k (b) - & k ( b') } } . (D5) 

And since 3§ k { b) = ±1, if \S§ k {b) + 38 k (h')\ = 2, then \3§ k {b) — £$ k (b')\ = 0, and vice versa. Consequently, using 
£^k(&) = ±1, it is easy to conclude that the absolute value of the above average cannot exceed 2, just as Bell concluded: 

- 2 < ^^4(a)^ fc (b) + £/ k (&) & k (b') + srf k (a!) @ k (b) - ,e4(a') ^fc(b') ^ < +2. (D6) 

Let us now try to understand why the replacement (ID4I) is illegal. To begin with, Einstein’s (or even Bell’s own) 
notion of local-realism does not demand this replacement. Since this notion is captured already in the definition of 
the measurement functions £/ (a, X k ), the LHS of (ID4I) satisfies the demand of local-realism perfectly well. To be sure, 
mathematically there is nothing wrong with a replacement of four separate averages with a single average. Every 
school child knows that the sum of averages is equal to the average of sums. But this rule of thumb is not valid in 
the above case, because (a, b), (a, b'), (a', b), and (a', b') are mutually exclusive pairs of measurement directions, 
corresponding to four incompatible experiments. Each pair can be used by Alice and Bob for a given experiment, for 
all runs 1 to n, but no two of the four pairs can be used by them simultaneously. This is because Alice and Bob do not 
have the ability to make measurements along counterfactually possible pairs of directions such as (a, b) and (a, b') 
simultaneously. Alice, for example, can make measurements along a or a', but not along a and a ' at the same time. 

But this inconvenient fact is rather devastating for Bell’s argument, because it means that his replacement (ID4jl is 
illegitimate. Consider a specific run of the EPRB experiment and the corresponding quantity being averaged in (ID4I) : 


#/ k (a)& k (b) + ^4(a) & k (b') + £/ k (a!) @ k (b) - £/ k (a!) gg k (b'). 


(D7) 
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Here the index k = 1 now represents a specific run of the experiment. But since Alice and Bob have only two particles 
at their disposal for each run, only one of the four terms of the above sum is physically meaningful. In other words, 
the above quantity is physically meaningless, because Alice, for example, cannot align her detector along a and a' at 
the same time. And likewise, Bob cannot align his detector along b and b' at the same time. What is more, this will 
be true for all possible runs of the experiment, or equivalently for all possible pairs of particles. Which implies that all 
of the quantities listed below, as they appear in the average (1D5|) . are un-observable and hence physically meaningless: 


£/i(a)&i(b) + £%i(b) — 3§i(h') , 

,e4(a)^ 2 (b) + ,s/ 2 (a) ^ 2 (b') + .^(a') J? 2 (b) - s^(a’) ^ 2 (b'), 
■s/ 3 (a)^ 3 (b) + ,s/ 3 (a)^ 3 (b') + £/ 3 (a')& 3 (b) - £f 3 (a') & 3 (b'), 
,fi4(a) ^(b) + ^(a) ^(b') + ^(a') ^(b) — ^(a') 3§±(yJ ), 


£f n (a)& n (b) + £/ n (a) 3§ n {b') + &f n {a!) 38 n {b) - &/ n (a') 3§ n {b’). 

But since each of the quantities above is physically meaningless, their average appearing on the RHS of (ID4I) . namely 

^ ^4(a) ^fc(b) + ^4(a) ^fc(b') + ^4(a')^fc(b) - M^a') ^k(b') j , (D8) 

is also physically meaningless. That is to say, no physical experiment can ever be performed - even in principle - that 
can meaningfully allow to measure or evaluate the above average, since none of the above list of quantities could have 
experimentally observable values. Therefore the innocuous looking replacement (El made by Bell is, in fact, illegal. 


On the other hand, it is crucial to note that 


each of the four averages appearing on the LHS of replacement (ID4I) , 


£(a, b) 
£(a, b') 
£(a', b) 


lim 

ra»l 


lim 

n» 1 


lim 

ra»l 




£/(a, X k ) &{b, \ k ) 
£/(a, X k ) &{b\ X k ) 
£?{a', X k ) S8{ b, X k ) 


^ ^4 (a) & k {h) ^ , 
4(a)^ fe (b') ^ , 
^■*4(a')<^fc(b) ^ , 


(D9) 

(DIO) 

(Dll) 


and £(a'. b') 


lim 

n> 1 


1 

- E ^( a '> A ") ^( b, > Afc ) 


^■«4( a!)3§ k {b') ^ , 


(D12) 


L fe = l J 

is a perfectly well defined and observable physical quantity. Therefore the bounds (1D3I) on their sum are quite harmless. 
These bounds of {—4, +4}, however, have never been violated in any experiment (indeed, nothing can violate them). 


In summary, Bell and his followers derive the upper bound of 2 on the CHSH string of averages by an illegal move. 
In the middle of their derivation they unjustifiably replace an observable, and hence physically meaningful quantity, 

^ ^4(a) 3§k{b) ^ + ^^4(a)^fc(b')^ + ^ £s/ k (a') & k (h) ^ £/ k (a') @ k {b') ^ , (D13) 

with an experimentally un-observable, and hence physically entirely meaningless quantity 

^^4(a)^ fc (b) + srf k (a) SS k {b') + &t k {a!) 3S k (b) - ^4(a') 3§ k {b') ^ . (D14) 

If they do not make this illegitimate replacement, then the upper bound on the CHSH string of averages is 4, not 2. 
It is mind-boggling why such a blatant mistake has been overlooked by the physics community for over 50 years 14]. 
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One may suspect that the above conclusion is perhaps an artifact of the discrete version, (|D2I) . of the expectation 
values £ (a, b). Perhaps it can be ameliorated if we considered the CHSH sum (ID 1 1) in the following continuous form: 

f &/(a, X)38(b, X)dp(X) + f s3(a, A)^(b', A) dp(X) + f £/(&', A) 38(h, X)dp(X) — f srf(af, X)38(W, X)dp(X), 
Ja J a J a J a 

(D15) 

where A is the space of all hidden variables A and p( A) is the probability distribution of A [2]. Written in this form, it 
is now easy to see that the above CHSH sum of expectation values is both mathematically and physically identical to 


£/(&, A) { 38{b, A) + 38{b', A) } + si (a', A) { 38(b, A) - 38{b', X) } dp{ A). 


(D16) 


ja L J 

But since the above two integral expressions are identical to each other, we can use the second expression without 
loss of generality to prove that the criterion of reality used by Bell is excessively restrictive compared to that of EPR. 


Proof : As we saw above, (|D16l) involves an integration over fictitious quantities like si(a, A) {38{b, A) ± 38(b', A)}. 
These quantities are not parts of the space of all possible measurement outcomes such as si (a, A), si (a', A), 38 (b, A), 
38(b', A), etc.] because that space — although evidently closed under multiplication — is not closed under addition. 
This is analogous to the fact that the set O := { 1, 2, 3,4, 5,6} of all possible outcomes of a die throw is not closed under 
addition. For example, 3 + 6 is not a part of the set O. But there is also a much more serious physical problem with 
Bell’s criterion of reality. As noted above, the quantities si (a, A) {38( b, A) ± 38 { b', A)} are not physically meaningful 
in any possible physical world, classical or quantum. That is because 38{ b, A) and 38 (W, X) can coexist with si (a, A) 
only counterfactually. If 38{ b, A) coexists with si{ a, A), then 38( b', A) cannot coexist with si( a, A), and vice versa. 
But in the proof of his theorem Bell presumes both 38{ b, A) and 38 { b', A) to coexist with si( a, A) simultaneously. 
That is analogous to being in New York and Miami at exactly the same time. But no reasonable criterion of reality 
can justify such an unphysical demand. The EPR criterion of reality most certainly does not demand any such thing. 


In conclusion, since the two integrands of (ID16I) are physically meaningless, the bounds of —2 and +2 on (ID15I) are 
also physically meaningless. They are mathematical curiosities, without any relevance for the question of local realism. 
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